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ABSTRACT 


Previous  work  has  shown  that  a  speed-independent  dynamic  fracture 
toughness  property  can  be  used  in  an  elastodynamic  analysis  to  describe 
crack  initiation  and  unstable  propagation  under  impact  loading.  In  this 
paper,  a  further  step  is  taken  by  extending  the  analysis  from  simple  labor¬ 
atory  test  specimens  to  treat  more  realistic  crack-structure  geometries.  A 
circular  cylinder  with  an  initial  part-through  wall  crack  subjected  to  an 
impulsive  loading  on  its  inner  surface  is  considered.  The  crack  is  in  a 
radial-axial  plane  and  has  its  length  in  the  axial  direction  long  enough  that 
a  state  of  plane  strain  exists  at  the  center  of  the  crack.  Crack  growth 
initiation  and  propagation  through  the  wall  is  then  calculated.  It  is  found 
that,  once  Initiated,  crack  propagation  will  continue  until  the  crack  pene¬ 
trates  the  wall.  Crack  arrest  within  the  wall  does  not  appear  to  be  possible 
under  the  conditions  considered  in  this  paper. 


NOMENCLATURE 

*  crack  length ,  an 

c^  Raylalgh  wav*  apaad,  a/sec 

c  bar  wave  speed,  m/sec 
o 

E  elastic  modulus,  Pa 

2 

G  dynamic  strain  energy  release  rate,  J/a 
h  cylinder  wall  thickness,  on 

1/2 

Kj  dynamic  stress  Intensity  factor,  MPam 

1/2 

Kja  static  crack  arrest  toughness,  MPam 

1/2 

Kid  dynamic  fracture  toughness,  MPsn 

1/2 

dynamic  initiation  toughness,  MPam 

1/2 

minimum  value  of  MPam 
R  cylinder  mean  radius,  mn 

T  period  of  oscillation,  u  sec 

t  time,  u  sec 

U  total  energy  per  unit  length,  J/mm 

V  crack  speed,  m/scc 


INTRODUCTION 

Fracture  mechanlca  offera  significant  opportuni¬ 
ties  for  developing  failure-safe  structures.  By  relat¬ 
ing  the  fracture-critical  flaw  sixes  at  various 
locations  in  the  structure  to  the  expected  applied 
stresses,  more  accurate  specification  of  both  material 
toughness  requlrsments  and  of  non-destructive  inspec¬ 
tion  limits  will  be  possible.  This  must  lead  towards 
an  effective  uae  of  materials  without  compromising  the 
Integrity  of  the  structure.  But,  first,  an  extension 
of  the  relatively  simple  fracture  mechanics  techniques 
that  are  now  available  must  be  made  to  treat  the  mater¬ 
ials,  geometries,  and  loading  conditlona  arising  in 
real  engineering  application. 

The  specific  objective  of  this  work  is  the  devel¬ 
opment  of  fracture  mechanics  analysis  procedures  for 
crack  growth  initiation,  unstable  propagation  and 
arrest  in  ship  structures  subjected  to  blast  loading. 
Because  there  are  several  different  aspects  of  this 
problem  that  are  beyond  current  analysis  capabilities, 
a  multi-step  approach  is  being  followed.  A  first  step 
was  to  assess  the  effect  of  impact  loading  in  a  com¬ 
bined  experimental  and  mathematical  analyals  program 
using  simple  laboratory  test  specimens.  The  work  re¬ 
ported  In  this  paper  extends  this  work  to  more  eomplex  ■ 
structures.  The  results  are  expected  to  permit  the 
effects  of  structural  geometry  on  dynamic  crack  propa¬ 
gation  to  be  assessed  and,  in  addition,  will  assist  in 
the  design  of  a  critical  experimental  test  of  the 
validity  of  dynamic  fracture  mechanics  predictions 
under  shock  loading  conditions. 

PRELIMINARY  DISCUSSION 

Dynamic  fracture  mechanics  encompasses  all  pro¬ 
blems  involving  crack  growth  initiation  and  subsequent 
unstsbls  propsgatlon  up  to  and  including  crack  arrest. 
The  methodology  was  developed  to  treat  problems  where, 
for  an  acceptable  solution,  inertia  forces  must  be 
Included  in  the  equations  of  motion  of  the  cracked 
body.  At  present,  dynamic  fracture  mechanics  treat¬ 
ments  are  limited  to  problems  where  the  basic 


assumptions  of  linear  alaatlc  fracture  mechanics 
(LEFM)  are  valid.  The  eaaential  assumption  in  this 
approach  is  chat  the  plastic  deformation  attending  the 
propagating  crack  tip  is  small  enough  to  be  ''domina¬ 
ted”  by  the  elastic  deformation  in  the  field  surround¬ 
ing  the  crack  tip.  In  these  conditions,  the  plastic 
energy  absorbed  in  the  fracture  process  can  be  taken 
as  characteristic  of  the  material  vlth  the  body  being 
treated  as  completely  elastic.  Problems  of  crack 
growth  initiation  and  subsequent  rapid  unstable  crack 
propagation  can  then  be  solved  by  using  elastodynaml- 
cally  computed  stress  intensity  factors  coupled  with 
experimentally  determined  dynamic  fracture  toughness 
values. 

The  stress  intensity  factor  Kj  enters  in  the  com¬ 
puted  elastodynamic  stress  field  in  the  immediate 
vicinity  of  the  crack  tip.  It  can  depend  on  time  t, 
the  crack  tip  speed,  V,  the  crack  length,  the  exter-  ; 
nal  geometry  of  the  cracked  body,  material  constants, 
and  on  the  applied  loads.  The  conditions  governing  j 
crack  motion  in  a  body  can  be  expressed  in  terms  or 
Kj(t,V)  and  experimentally  determined  critical  values 
that  are  taken  to  be  properties  of  the  material. 

Thus,  for  a  propagating  crack 

Kx(t,V)  -  KjjjOO,  (1) 

where  K._  is  known  as  the  dynamic  fracture  toughness. 
Ordinarily,  KT_  values  will  be  greater  than  al¬ 
though  it  is  possible  that  K^,  the  minimum  value  of 
K_  ,  can  be  less  than  K„. 

An  equality  is  sometimes  used  for  crack  arrest. 
This  is  expressed  in  terms  of  Kj  and  an  "arrest  tough¬ 
ness"  parameter,  K.  .  However,  while  the  concept  can 
be  useful  as  an  appfoxlmatlon,  crack  arrest  is  more 
rigorously  defined  as  occurring  only  when  Equation  (1) 
cannot  be  aatiafied.  That  is,  the  crack  will  arrest 
at  a  time  t  when  R-  <  K-  for  all  t  >  t  .  Thus, 
crack  arrest  is  properly  viewed  as  the  termination  of 
a  general  dynamic  crack  propagation  process,  not  aa  a 
unique  event  governed  solely  by  material  properties. 
This  is  the  way  in  which  the  methodology  will  be 
applied  here. 

Because  of  the  equivalence  that  exists  between 
them,  applications  of  LEFM  can  be  made  either  in  terms 
of  the  stress  intensity  factor  or  the  energy  release 
rate  parameter  6.  That  is,  for  plana  strain  condi¬ 
tions 

2 

G  -  ±2-  A(V)  ,  (2) 


ture  toughness  in  which  a  quasi-static  calculation 
would  predict  crack  arrest  at  a  penetration  equal  to 
71Z  of  the  wall  thickness,  the  dynamic  fracture  me¬ 
chanics  calculation  of  Reference  3  revealed  that  the 
crack  would  completely  penetrate  the  wall. 

The  growth  of  a  part-through  wall  crack  is  gener¬ 
ally  a  three-dimensional  problem.  An  effective  sim¬ 
plifying  approach  is  to  assume  a  plane  strain  ideali¬ 
zation;  e.g.,  considering  an  initial  part-through  wall 
flaw  to  have  a  much  greater  length  in  the  direction 
parallel  to  the  surface  of  the  wall  than  in  the  thick¬ 
ness  direction.  This  approach  was  successfully  taken 
in  the  finite-difference  computation  for  a  run-arrest 
event  in  a  thermally  shocked  pressure  vessel  with  an 
internal  part-through  wall  crack  renorted  in-  Reference 
(3).  This  method  has  been  extended  here  to  take 
account  of  rapidly  applied  loading  for  application  to 
structures  under  impact  loading. 

The  following  analyses  have  two  primary  purposes. 
The  work  reported  in  Reference  (1)  has  demonstrated 
that  the  elastodynamic  crack  propagation  and  crack 
arrest  methodology  can  be  applied  to  impact-loaded 
dynamic  tear  test  specimens.  It  is  therefore  desir¬ 
able  to  study  rapid  fracture  and  crack  arrest  in  an 
impulsively  loaded  component  that  is  more  realistic 
than  a  simple  laboratory  specimen.  Second,  because 
there  is  a  lack  of  Information  on  atresa  intensity 
factors  for  impulsively  loaded  structures,  these 
analyses  will  provide  useful  quantitative  predictions 
of  the  magnitude  of  the  impulsive  loading  required  to 
initiate  crack  propagation  in  a  precracked  or  flawed 
cylindrical  vessel.  It  is  expected  that  this  informa¬ 
tion  can  be  used  to  design  a  critical  test  of  the 
methodology. 

ANALYSIS  AND  RESULTS 

The  Computational  Model 

The  objectives  of  immediate  Interest  are  to  deter¬ 
mine  the  magnitude  of  the  impulsive  loading  that  will 
Initiate  propagation  of  a  part-through  crack  in  the 
wall  of  a  cylindrical  vessel  and  to  determine  whether 
or  not  this  crack  will  arrest  before  it  penetrates  the 
wall.  The  particular  loading  envisioned  here  is  an 
Intense  pressure  spike.  The  duration  of  the  spike 
will  be  such  smaller  than  the  fundamental  period  of 
oscillation  of  the  cylinder  to  constitute  an  Impulsive 
or  a  shock  loading. 

The  cylindrical  vessel  to  be  addressed  in  the 
following  is  depicted  in  Figure  1.  It  consists  of  a 


where  A(V)  is  a  universal  geometry-independent  func¬ 
tion  that  is  unity  at  zero  crack  speed,  and  Increases 
monotonleally  to  become  unbounded  mV*  c_,  where  c. 
is  the  speed  of  Rayleigh  waves  in  the  material.  For 
most  practical  altuatlona,  A(V)  can  be  taken  as  equal 
to  one. 

In  addition  to  the  work  reported  by  Kanninan,  et 
al  (1),  other  applications  which  illustrate  the 
necessity  of  a  dynamic  approach  have  been  given  by 
Hahn,  et  al  (2.3).  Specifically,  a  series  of  thermal 
shock  experloentt  performed  on  axially-cracked  thlck- 
walleJ  cylinders  at  the  Oak  Ridge  National  Laboratory 
(4)  were  analyzed.  It  was  found  that  the  predictions 
were  in  quite  good  agreement  vith  the  measured  dis¬ 
tances  of  crack  propagation  prior  to-erresg.  These 
experiments  produced  only  short  crack  jump  lengths 
(e.g.,  penetration  to  14Z  of  the  wall  thickness)  and 
for  these  a  quasi-static  approach  would  also  be  accep¬ 
table.  However,  for  conditions  which  would  produce 
much  larger  penetrations,  a  significant  difference  was 
found.  Specifically,  for  a  postulated  reduced  frac- 


Flgure  1.  Cracked  Cylinder  with  Finite  Difference  Grid 


vary  long  circular  cylladar  of  mu  radius  R  and  wall 
thickness  h.  A  part-through  wall  crack  of  dapth  aQ 
la  aaauaad  to  azlat  at  ths  inalda  aurfaca  of  the 
cylinder.  The  length  of  the  crack,  like  Che  cylinder, 
la  taken  to  be  long  enough  for  plane  strain  conditions 
to  be  valid.  This  type  of  crack  clearly  poses  a  more 
severe  condition  than  one  with  a  finite  axial  lsngth. 
An  internal  flow  is  considered. 

The  loading  is  a  unifora  internal  iapulslve  pres¬ 
sure.  The  reason  for  selecting  an  internal  pressure 
over  an  external  pressure  Is  also  one  of  convenience. 
The  Internal  loading  iaaediately  causes  the  crack  to 
open.  By  contrast,  if  an  external  lapulsa  were 
applied,  the  cylinder  would  have  to  experience  half  an 
oscillation  before  any  crack  opening  would  occur. 
Hence,  ths  Internal  loading  requires  less  computer 
tine.  Also,  the  danger  of  buckling  the  cylinder  la 
reduced. 

This  analysis  is  restricted  to  conditions  where 
the  time  required  for  a  longitudinal  wave  to  propagate 
through  the  wall  Is  small  compared  to  the  fundamental 
period  of  oscillation  of  the  cylinder.  Hence,  it 
applies  when  the  thickness  of  the  cylinder  is  small 
compared  to  Its  radius.  However,  even  though  the 
cylinder  Is  a  thin  shell,  it  does  not  follow  that 
classical  shell  theory  Is  applicable  to  the  present 
problem. 

The  net  effect  of  an  Impulsive  loading  on  a  thin 
cylinder  la  to  Impart  s  nearly  uniform  radial  velocity 
v  to  Its  walls.  (This  condition  is  frequently 
aSsumed  when  analyzing  Impulsive  loading  of  thin 
shells.)  Under  such  a  loading,  the  cylinder  will 
respond  by  oscillating  in  essentially  the  breathing 
mode  with  superimposed  low  frequency  flexural  modes. 
Had  an  initial  velocity  Imparted  only  to  the  Inside 
surface  been  considered,  then  higher  frequency  modes 
would  have  been  excited.  Computations  for  the  latter 
case  showed  that  the  period  of  these  higher  modes  Is 
small  compared  to  the  fracture  event.  While  they  can 
be  Included,  they  do  little  more  then  cloud  the  under¬ 
standing  of  the  fracture  phenomenon.  Furthermore, 
these  high  frequency  modes  would  be  the  first  to  be 
damped  out  and,  hence,  can  also  be  neglected  on  this 
basis. 


of  motion  are  solved  with  stepwise  increasing  time 
with  the  fracture  criterion.  Equation  (1),  being  tested 
in  each  time  step.  When  the  fracture  criterion  la 
satisfied,  the  crack  la  permitted  to  advance  one  half 
of  the  radial  nodal  spacing.  For  the  computations 
described  herein,  the  number  of  nodes  through  the 
thickness  was  kept  constant  at  thirteen.  The  number 
of  nodes  in  the  circumferential  direction  was  selected 
to  keep  the  aspect  ratio  (the  ratio  of  the  radial  grid 
dimension  to  the  circumferential  grid  dimension)  of 
the  grid  at  approximately  0.04.  Favorable  results 
have  been  obtained  previously  for  this  aspect  ratio. 

Computational  Results  for  a  Stationary  Crack 


Figure  2  shows  the  variation  of  the  calculated 
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Figure  2.  Stress  Intensity  Factor  Versus  Time  For  A 
Stationary  Crack 

stress  intensity  factor  with  time  for  an  impulsively 
loaded  steal  cylinder  for  a  stationary  crack.  These 
results  are  for  R  -  305-mm,  h  -  26-mm,  a  -  6-mm,  and 
v  -  lm/s.  Note  that  the  Intervals  where  K  •  0  cor¬ 
respond  to  the  closing  of  the  crack.  Of  most  impor¬ 
tance  Is  the  maxlsma  stress  intensity  factor  - 

9.5  MNm~3/2  in  this  instance)  determined  In  this  way. 
This  value  is  required  in  order  to  determine  the  mini¬ 
mum  load  necessary  to  initiate  crack  propagation. 

The  period  T  of  oscillation  in  the  breathing  mode 
of  an  unflawed  cylinder  Is 

,  _  2*R(1  -  v2)1/2 


The  Solution  Frocedure 

In  common  with  the  work  presented  in  Reference 
(1),  the  approach  la  within  the  confines  of  LEFM  with 
Inertia  effects  Included.  The  finite  difference 
method  Is  used  to  integrate  the  equations  of  notion 
expressed  In  terms  of  the  nodal  displacements2.  The 
finite  difference  grid  is  depicted  in  Figure  1.  It 
can  be  seen  that,  because  of  the  symmetry  with  respect 
to  the  diametrical  plana  containing  tha  crack,  only 
half  the  cylinder  need  be  analyzed.  The  faces  of  the 
crack  nay  opan  up  but  they  are  prevented  from  pene¬ 
trating  each  other.  The  stress  intensity  factor  la 
determined  from  (2)  using  the  energy  release  rate  cal¬ 
culated  from  the  displacements  at  the  nodes  In  the 
near  vicinity  of  the  crack  tip. 

At  tlae  t  •  0,  each  node  is  given  only  a  radial 
initial  velocity  v  .  The  Impulse  per  unit  of  surface 
area  la 

I  ■  ohv  (3) 

o 

where  p  is  ths  density  of  the  material.  The  equations 


where  C  -  /e/p  is  the  bar  wave  speed  which  is  5000 
m/s  for  steel.  For  v  »  0.3,  the  period  for  R  -  305-nm 
Is  365us.  This  Is  very  nearly  equal  to  tha  period  of 
the  stress  intensity  factor  In  Figure  2.  It  Is  also 
apparent  from  this  figure  that,  while  the  amplitude 
varies  slightly.  It  can  be  taken  as  essentially  con¬ 
stant.  Therefore,  subsequent  analyses  of  the  station¬ 
ary  crack  need  only  consider  the  first  half  period  of 
oscillation  to  determine  reasonable  values  of  K  . 

The  maximum  stress  in  an  unflawed  cylinder  Sub¬ 
jected  to  a  uniform  Impulsive  loading  is  readily  shown 
to  be 


C0(l  -  v2)1/2 


which  is  Independent  of  R/h.  The  maximum  stress  in¬ 
tensity  factor  can  be  related  to  o  by  a  relation 
valid  In  quasi-static  conditions.  °Thls  is 


/sa“  f(a  /h,  R/h) 
o  0 


2  A  detailed  description  of  the  method  appears  in  Raf- 
erencs  (2)  with  an  abbreviated  account  given  In  (3). 


where  f  is  a  function  of  the  geometry  of  the  flawed 
cylinder.  In  Table  1  the  computed  K  values  for 
various  geometries  are  shown.  Thesenlve  also  been 


Table  1.  Computational  Results  For  A  Stationary  Crack 
In  An  Impulsively  Loaded  Radially  Cracked 
Cylinder 
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11.7 
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1.78 

1.54 

30. 

14.5 

62.7 

1.85 

1.55 

normalized  with  respect  to  c  <Va  .  It  appears  from 
these  results  that  the  normalized  stress  intensity 
factor  is  virtually  Independent  of  R/h  for  these  thin 
walled  cylinders.  Hence,  for  design  purposes,  it 
could  be  considered  that 


<l-v2)1/2 


o  .  .1/2 

Co  0 


where  the  result  given  in  Table  1  has  been  used. 

The  maximum  stress  intensity  factor  also  appears 
to  be  nearly  independent  of  a  /h.  However,  only  rela¬ 
tively  small  values  of  a  /h  were  used  here  and  it  is 
likely  that  the  edge  effect  is  small  in  this  range. 
Results  for  steady  state  vibrations  of  center-crack 
plates  and  an  infinite  plane  with  a  periodic  system 
of  cracks— see  Reference  (S)— also  indicate  that  edge 
effects  do  noc  become  significant  until  ag/h  >  0.6. 

For  larger  values  of  a  /h,  IC^/cj  /ira_  would  be  expec¬ 
ted  co  be  greater  than°the  value  found  hare. 

Table  1  also  shows  the  maximum  stress  attained  at 
the  point  two  nodes  (4-mm)  ahead  of  the  crack  tip. 

This  value  is  denoted  as  amax-  Hare  agsln,  there  is 
only  a  slight  dependence  on* R/h.  But,  o  does  de¬ 
pend  upon  a  /h.  By  comparing  o  to  the  yield 
stress,  these  results  can  be  used* Co  predict  Che  maxi¬ 
mum  Impulse  chat  would  satisfy  the  conditions  of  LEFM. 

Computational  Results  for  Unstable  Crack  Propagation 

For  the  same  Initial  loading,  the  computations 
for  a  stationary  crack  shov  that  K  increases  with 
increasing  crack  depth.  If  crack  propagation  were  to 
initiate,  it  would  therefore  appear  that  the  crsck 
would  be  propagating  in  an  Increasing  K  field.  Under 
such  circumstances  the  crack  would  not  arrest,  but 
would  certainly  penetrate  the  wall.  On  the  other 
hand,  the  impulsive  loading  only  Imparts  a  finite 
amount  of  energy  co  the  cylinder.  And,  as  the  crack 
propagates,  it  must  consume  some  of  this  energy  in  the 
fracture  process.  Hence,  there  is  a  question  as  to 
whether  or  not  this  energy  loss  is  sufficient  to  re¬ 
duce  the  amplitude  of  the  stress  intensity  factor  suf¬ 
ficiently  to  arrest  the  crack.  This  question  can  be 
properly  addressed  only  with  a  dynamic  analysis. 

A  quarter-through-wall  crack  is  frequently  taken 
to  be  minimum  identifiable  crack  length.  If,  for 


reasons  of  safety,  cracks  exceeding  half  the  wall 
thickness  are  deemed  not  permissible,  a  question  of 
concern  is  whether  or  not  a  quarter-through  wall  crack 
which  initiates  will  arrest  before  it  propagates  half 
way  through  the  wall.  To  study  this  question,  assume 
that  the  material  has  a  speed  Independent  toughness 
kID  "  Kic  ud  the  amplitude  of  the  final  stress  inten¬ 
sity  factor  is  K  •  i.e.,  the  crack  just  arrests.  A 
balance  of  energy  per  unit  length  requires  that 

k|„(1  -  v2) 

»»  11  _  _ /  _  _  V  /AV 


where  U  denotes  the  total  energy  per  unit  length  and 
the  subscripts  o  and  f  are  used  to  denote  initial  and 
final  quantities.  Using  (2),  the  energies  can  be 
written  as 
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irpv  Rh 
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^(1  -  S)  Rh 
3.2  E  a 


K2dU  -  v2)  Rh 
3.2  E  a. 


5  Ti  1/2  • 


Substituting  Equations  (9)  into  (8)  with  a  /h  ■  1/A  and 
a./h  “  1/2  gives  ° 


-  (0.20  j  ♦  .5)1/2 


But,  for  the  crack  to  initiate,  Rid  <  (K  )  .  How¬ 

ever,  this  is  inconsistent  with  (10)and  the^ypothesls 
that  the  crack  arrests.  Therefore,  this  contradiction 
Implies  that  the  crack  will  not  arrest  before  it  propa¬ 
gates  half-way  through  the  wall.  Because  of  the  1m- 
pect  loading,  the  crack  might  initiate  with  K- .  £ 

(Kq)  x  <  Kj_  and  no  Inconsistency  would  appear.  Also, 
asnoced  earlier,  viscous  damping,  .which  could  have  a 
significant  effect  under  certain  circumstances,  is  not 
Included. 

In  order  to  determine  the  character  of  crsck  prop¬ 
agation  in  a  cylinder  under  impulsive  or  shock  loading, 
computations  were  performed  for  a  steel  cylinder  of 
radius  R  -  305-mm  and  thickness  h  ■  26-mm  having  an 
initial  radial  crack  of  depth  a,,  -  6 -mm  (a^h  -  0.23). 
An  initial  uniform  radial  velocity  of  lm/s  was  imparted 
to  the  cylinder.  These  are  the  same  conditions  upon 
which  the  results  of  Figure  2  are  beaed.  The  fracture 
toughness  was  taken  as  a  speed  independent  value  equal 
to  98X  of  the  maximum  stress  intensity  factor  experi¬ 
enced  by  a  stationary  crack  for  chase  sane  conditions; 
i.e.,  KiD  ”  9.25  MN/m3/2  per  unit  of  initial  velocity. 

A  plot  of  computed  crack  length  versus  time  is 
shown  in  Figure  3.  It  can  be  seen  that  approximately 
84  u  sec  are  required  for  the  stress  intensity  to 
build  up  to  the  critical  value  and  to  initiate  crack 
growth.  During  this  Interval,  the  variation  of  the 
stress  Intensity  factor  with  time  is  depicted  in  Fig¬ 
ure  2.  After  initiation,  the  crack  propagates  at  a 
speed  of  approximately  300  m/s.  As  the  crack  tip 
propagated  further  into  the  wall,  its  speed  increasea. 
This  reflects  the  unstable  nature  of  the  crack  gTowth. 
Some  20  u  sec  after  initiation  the  crack  tip  finally 
penetrated  the  exterior  surface.  The  speed  of  the 
crack  tip  during  the  final  stages  of  the  event  was 
approximately  2200  m/s. 


Figure  3.  Crack  Propagation  In  Impulsively  Loaded 
Crack  Cylinder 


For  all  practical  purpoaae  chla  computation  con¬ 
fine  the  earlier  auaplcion  that,  if  growth  of  an 
internal  crack  commences,  the  crack  tip  will  eventu¬ 
ally  penetrate  the  exterior  surface.  This  conclusion 
is  predicated  upon  the  assumption  that  the  material  la 
elastic-brittle  and  has  a  crack-speed-independent 
fracture  toughness.  In  reality,  as  the  remaining 
ligament  becomes  sufficiently  small,  inelastic  be¬ 
havior  becomes  important.  Then,  the  fracture  process 
is  no  longer  K-domlnated  as  required  for  LEFM  to  be 
applicable.  While  such  inelastic  behavior  is  impor¬ 
tant  and  of  Interest,  its  considerstion  is  beyond  the 
scope  of  this  paper  and  is  reserved  for  further 
research. 

CONCLUSIONS 

Elastodynamlc  crack  propagation  calculations  have 
been  performed  for  a  part-through  cracked  circular 
cylinder  subjected  to  a  uniform  impulsive  losding. 
These  calculations  have  demonstrated  that,  while  ac¬ 
tual  crack-structure  geometries  are  more  complicated 
than  are  the  simple  laboratory  test  specimens  pre¬ 
viously  analyzed,  they  can  be  treated  effectively  with 
dynamic  fracture  mechanics.  Specifically,  for  the 
particular  geometry  considered  in  this  paper: 

1.  The  minimum  impulse  required  to  initiate 
unstable  crack  propagation  for  a  given 
crack  depth  has  been  determined  in  terms 
of  material,  mechanical,  and  fracture 
properties. 

2.  Once  initiated,  crack  propagation  in  the 
through-wall  direction  will  continue 
until  the  crack  penetrates  the  wall- 
crack  arrest  within  the  wall  does  not 
appear  to  be  possible. 

3.  Crack  length  time  predictions  can  be 
made  for  comparison  with  experimentally 
measured  times  of  crack  growth  initia¬ 
tion  and  of  crack  penetration. 


should  be  recognized  that  the  specific  conclusions 
cited  are  restricted  to  LEFM  conditions  and  for  a 
material  with  a  relatively  conatant  dynamic  fracture 
toughness. 
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It  is  clear  that  more  complicated  geometries  can  be 
addressed  with  this  anslysis  procedure  (vlth  a  corres¬ 
ponding  Increase  la  computer  costs).  However,  it 


